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An extreme ultraviolet (EUV) single attosecond pulse passing through a laser-dressed dense gas 
is studied theoretically. The weak EUV pulse pumps the helium gas from the ground state to the 
2s2p( 1 P) autoionizing state, which is coupled to the 2s 2 ( 1 S) autoionizing state by a femtosecond 
infrared laser with the intensity in the order of 10 12 W/cm 2 . The simulation shows how the transient 
absorption and emission of the EUV are modified by the coupling laser. A simple analytical ex- 
pression for the atomic response derived for 5-function pulses reveals the strong modification of the 
Fano lineshape in the spectra, where these features are quite universal and remain valid for realistic 
pulse conditions. We further account for the propagation of pulses in the medium and show that the 
EUV signal at the atomic resonance can be enhanced in the gaseous medium by more than 50% for 
specifically adjusted laser parameters, and that this enhancement persists as the EUV propagates 
in the gaseous medium. Our result demonstrates the high-level control of nonlinear optical effects 
that are achievable with attosecond pulses. 

PACS numbers: 32.80.Qk, 32.80.Zb, 42.50.Gy 



I. INTRODUCTION 

The rapidly developing technologies of ultrashort light 
sources have gained wide applications in the past 10 
years [l], most notably with the use of attosecond 
pulses in time-resolved spectroscopy. A typical measure- 
ment scheme shines an extreme ultraviolet (EUV) light 
pulse, in the form of an attosecond pulse train (APT) or 
a single attosecond pulse (SAP), together with a synchro- 
nized infrared (IR) laser pulse, on an atomic or molecular 
target, and measures the photoions, photoelectrons, or 
photoabsorption of the EUV. The adjustable time delay 
between the EUV and the IR gives information of the 
dynamics of the target. For the SAP in the x-ray regime, 
this technique has been used to study the emission of 
Auger electrons Q in the so-called "streaking" model. 
Down to the EUV energy, a similar scheme was used to 
time-resolve autoionization where the IR primarily just 
depleted the resonances excited by the EUV 4, 5]. For a 
SAP at even lower energies, a series of bound states can 
be excited at once. These states can then be ionized by 
the IR through different quantum pathways to exhibit in- 
terference patterns , or dressed by the IR temporarily 
to exhibit the ac Stark shift Q- 

The studies of coupled autoionizing states (AIS) by 
lasers have been carried out over the past 30 years theo- 
retically and experimentally [12|. In these earlier 
investigations, long pulses were used. For a typical AIS 
with a lifetime up to tenths of femtoseconds, the coupling 
lights were considered to be monochromatic, and the 
measured spectrum was obtained by scanning the photon 
energy over the widths of the resonances. Not until a few 
years ago was the same coupling scheme extended to the 
EUV energy range with femtosecond pulses, where tran- 
sient absorptions were measured or calculated [l3T - [l5j ]. 



Very recently, with the emerging attosecond pulses, the 
timcscalc of the coupling has been pushed shorter than 
the decay lifetime of an AIS d, [|| for direct observation 
in the time domain. In the same scheme, numerous new 
spectroscopic features were suggested by simulations [l6r - 
[211 ] . Some of these features have been shown recently by 
experiments with improved energy resolution |22| of the 
spectrometers. Further technological developments will 
likely bring more discoveries. 

In this work, we focus on resonant laser coupling be- 
tween two AISs, where an SAP passing through the 
medium is strongly reshaped in its spectral and temporal 
distributions by controlling the parameters of the laser 
and the medium. This reshaping of the SAP, which be- 
haves much beyond what can be described by the simple 
absorption rate or Beer's law [23|, has not been investi- 
gated so far. In Ref. [3], the actual propagation of EUV 
pulses was studied; however, those pulses with durations 
of tens of femtoseconds are too narrow in bandwidth to 
demonstrate any meaningful modification in the resonant 
spectrum. In this study, we look for how an SAP evolves 
in the medium, which cannot be revealed by a single- 
atom absorption cross section. 

For demonstration at the single-atom model, a 200 as 
EUV pulse and a time-delayed 9-fs laser pulse are applied 
to the helium atom. The 2s2p( 1 P) and 2s 2 ( 1 S) AISs 
are coupled by a 540-nm laser with intensities between 
10 12 and 10 13 W/cm 2 , where the weak EUV excites 2s2p 
from the ground state. For clarity, the pulses are not 
distinguished as the pump or the probe since they just 
couple different sets of states. In the single-atom calcu- 
lation, the electronic wave function, photoelectron, and 
EUV absorption spectra are carried out for a given set 
of pulses and atomic parameters. In order to extract the 
universal features resulting from the coupling, we further 



2 



derive simple analytic forms of the spectra for short light 
pulses that are approximated by 6 functions in time. This 
is a valid approximation if the pulse durations are shorter 
than the atomic timescales, e.g., the 17 fs decay lifetime 
of 2s2p. The Fano q parameter [2J| preserved in the final 
forms of the spectra indicates how the population trans- 
fer between the AISs is controlled by the coupling and 
how the wave packet evolves. 

For the macroscopic model, the pulses are allowed to 
propagate through a 2 mm helium gas with number den- 
sity 8 x 10 -3 cm -3 (or a pressure of 25 Torr at room tem- 
perature). The detected EUV spectra as the pulses exit 
the medium are calculated. We found that at certain cou- 
pling and gas conditions, an enhancement of more than 
50% of the incident light at the resonance is produced in 
the transmission spectrum. This enhancement peak in 
the frequency domain persists through the propagation 
where the surrounding background decays along the light 
path. This demonstrates a remarkable extension of the 
nonlinear optical control to the attosecond time regime. 
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FIG. 1: (Color online) Coupling scheme used in the model. 
The thick red arrow is the laser coupling. The thin purple 
arrow represents the EUV transition from the ground state 
to the resonance, where the bandwidth of the EUV 

pulse, indicated by the vertical purple line, widely covers the 
whole lineshape of the resonance. The green arrows with hol- 
low heads are the configuration interactions responsible for 
the autoionization. 



II. MODEL 

A. Single-atom wave function 

In this section, we derive the time-dependent total 
wave function of a three-level autoionizing system cou- 
pled by an EUV pulse and a laser pulse. Two AISs, 
composed by the bound parts \bi) and I&2), and the asso- 
ciated background continua \E\) and \E 2 ), respectively, 
are coupled by the laser, while the |6i)-|i?i) resonance 
is coupled to the ground state \g) by the EUV, both 
by dipole transitions. The coupling scheme can be ei- 
ther 3 type or A type. A schematic of the system is 
in Fig. [T] The pulses are linearly polarized in the same 
direction and collinearly propagated. The EUV field is 
expressed by E x (t) = F x (t)e lulxt + F x {t)e- lUxt where 
Fx(t) is the envelope and ujx is the central frequency 
of the pulse. The pulse envelope is generally a complex 
function, which includes all the phase factors in addition 
to the carrier- frequency terms. The electric field E^t) of 
the laser pulse, however, is kept in the exact form since 
it may be a few-cycle pulse. The total wave function is 
in the general form of 



m)) = e-^c g {t)\9) 



-C b2 (t)\b 2 



C&1WIM+ / c El {t)\Ei)dE 1 



c E2 (t)\E 2 )dE 2 



(1) 



In solving the time-dependent Schrodinger equation 
(TDSE), we make the following approximations: (1) 
All the second-order (two-electron) dipole matrix ele- 
ments, except the resonant excitation (bi\fi\g) (e.g., the 
transition between 2s2p and Is 2 ), are neglected, i.e., 



(.E2H&1) = (&2H-E1) = 0. (2) The rotating wave ap- 
proximation (RWA) is applied for the EUV coupling 
(but not the laser coupling) for its high frequency and 
low intensity. (3) The matrix elements of the Hamilto- 
nian involving \E\) and I-E2) are independent of energy, 
where we are only concerned with the energy ranges near 
the resonances. This means that the target structure is 
given in terms of the resonance energies E bl and Ef> 2 , the 
widths Ti and T2, and the q parameter of the 
resonance, which together constitute the basic Fano line- 
shapes [241 ] . Note that, with our first approximation, q 
cannot be defined for the 1 62 )- 1 ^2 ) resonance. 

Within the prescribed truncated space, the TDSE for 
the total wave function in Eq. (fT]) is reduced to the cou- 
pled equations of the coefficients therein, where these co- 
efficients as functions in time should be calculated ex- 
actly with the given initial conditions. However, the co- 
efficients associated with the continua, E\ and E 2 , make 
the calculation less tractable. Thus, we adopt the adia- 
batic elimination of the continua by assuming that they 
change much more slowly than the bound states. Then, 
the approximation CE 1 (t) = CE 2 {t) = in the coupled 
equations gives 



c El (t) 



Si \ Ei [-Fx(t)D* gE Mt) + ViC bl (t)} (2) 

c E 2 (t) = — —^— V2Cb 2 (t), 



(3) 



where 5i = Ex — E bl and 8 2 = Ex — Eb 2 are the detuning 
of the EUV relative to the transitions from the ground 
state to the AISs, AEi =E X -E a and AE 2 = E 2 ~E b are 
the continuous energies relative to the resonance energies, 
the Ds are dipole matrix elements, and V\ = {Ex\H\b\) 
and V 2 = (E 2 \H\b 2 ) are the strengths of the configura- 
tion interaction. In return, the coupled equations for the 
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bound state coefficients are reduced to satisfy 

ic g (t) = - F x (t)D* gbi Xc bl (t) - i \F x (t)\ 2 j gg c g (t) (4) 
ic bl (t) = - F x (t)D gbl \c g (t) 

- (Si +ini)c bl (t) - E L (t)D blb2 c b2 (t) (5) 

ic b2 (t) =- (S 2 + in 2 )c b2 (t) - E L (t)Dl lb2 c bl (t), (6) 

where k,\ = Ti/2 and k 2 = T 2 /2 are half widths of the 
resonances, A = 1 — i/q represents the transition from the 
ground state to \bi) and \E±) jointly, and j gg = ^DgE^ 2 
is proportional to the laser broadening. The dipoles (D) 
in this work are all real numbers since real orbitals are 
used for the bound states and standing waves are used for 
the continua. With Eqs. ((I])-®, the bound state part of 
the wave function is carried out. The procedure so far is 
identical to what has been employed in the early studies 
concerning coupled AISs |0l— fllj| , except that we do not 
apply RWA on the laser pulse. 

For pulses that arc long compared to the resonance life- 
time, the pulse bandwidth is narrow, thus each measure- 
ment gives the ionization yield at just one single photon 
energy in the spectrum. In this case, calculation of the 
ground-state population c g (t) described above would be 
enough for carrying out the ionization spectra. However, 
in the present study, the bandwidth of the attosecond 
pulse is broad and can excite the whole AIS, including 
a significant fraction of the background continuum. The 
continuum part of the wave function is needed in de- 
termining the electron or absorption spectra in a single 
measurement. Here, the coefficients of the continuum 
part are retrieved from the original coupled equations 

ic El {t) =(Ei - E x )c El (t) + V lCbl (t) 

iceM =(E 2 - E x )c E2 (t) + V 2 c b2 {t). (8) 



1. Photoelectrons 

In standard scattering theory, the atomic scattering 
waves for momentum k are given by 

= y^£^^^>W^:„(fc) (9) 

Irn 

in the energy-normalized form, where ui are real stand- 
ing wave radial functions, and rji are the scattering phase 
shifts. Both ui and r\i are determined by the atomic po- 
tential. If one detects photoelectrons, the momentum 
distribution will be the projection of the total wave func- 
tion at a large time onto the scattering waves, i.e., 

P(k)=Vm \{^(t))\\ (10) 

Likewise, if one measures the energy of the pho- 
toelectrons, the spectrum is given by P(E) = 
limt-j.oo |(?/>£ | \I/(i))| 2 , where are the scattering waves 
for energy E. Remember that \E\) and I-E2) are conven- 
tionally chosen as standing waves, which contains both 
the incoming and the outgoing components. At large 
times, since physically the AISs will decay to the con- 
tinuum, the wave packet will not have incoming com- 
ponents, i.e., the projection onto scattering waves is 
the same as the projection onto standing waves. Thus, 
the photoelectron spectra corresponding to the two res- 
onances are simply 

P x {Ei) = lim \c El (t)\\ 

t— voo 

P 2 (E 2 ) = lim \c E2 (t)\\ (11) 
t— ^ 00 

which require us to calculate the total wave function for a 
physical time much longer than the decay lifetimes, until 
Cs 1 (t) and CE 2 (t) stop evolving other than an oscillating 
phase. 



The retrieval of the continua can be viewed as a "cor- 
rection" after an iteration from the preliminary forms in 
Eqs. ([2]) and <j3j> - Now, the time-dependent total wave 
function in the form of Eq. ([T} is completely solved. 



B. Single-atom spectrum 



In the following, the electron energy spectrum and pho- 
toabsorption spectrum for a single atom are defined by 
the probability density per unit energy for emitting an 
electron and for absorbing a photon, respectively. In 
macroscopic cases where propagation is taken into ac- 
count, the spectroscopy is expressed in terms of the trans- 
mitted light intensity versus the photon energy. 



2. Photoab sorption 

In a given external field, the total energy absorbed by 
the atom is 



AU 



uiS(uj)duj, 



(12) 



where the response function S(uj) for ui > represents 
the absorption probability density, and is proportional to 
the pulse intensity in the frequency domain (T3 |. Thus, 
P(E) and S(co) have the same dimension and can be 
directly compared. The absorption cross section <j(uj) is 
related to S(uj) by 



'iiraujS^u)) 



E{lu) 



(13) 
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where a is the fine structure constant, and E(ui) is the 
Fourier transform of the electric field. By considering 
dipole interaction between the electric field and the atom, 
S(oj) is given by 



S(u) = -21m p,(uj)E*( 



(14) 



where /x(w) is the Fourier transform of the dipole mo- 
ment. 

With the wave function in Eq. ((1) , the dipole moment 
is given by 



(15) 



where the carrier frequency of the EUV has been factored 
out, and 

u x (t) = D gbl \*c* bl (t)c g (t) - iF x {t)j gg \c g (t)\ 2 , (16) 
u L (t) = D blb2 c* b2 (t)c bl (t). (17) 

Note that this model aims to deal with both short and 
long pulse lasers. In the case of a few-cycle laser, the 
Fourier transform p(ui) at very low frequencies may con- 
tain the contributions from both «l(£) and its complex 
conjugate, where the separation of the carrier frequency 
and the envelope is not beneficial; we thus do not employ 
this separation for the laser pulse in our study. 

In a typical ultrafast EUV-plus-IR experiment, the 
electron or absorption spectra measured over variable 
time delays between the two pulses are the most com- 
mon. Photoions are also usually measured, but they do 
not carry additional information. To study dynamic be- 
havior of autoionization directly in the time domain, light 
pulses shorter than the typical decay lifetime are neces- 
sary. By now few-cycle lasers of durations less than 10 fs 
are quite common. In the case where both pulses are sig- 
nificantly shorter than any other atomic timescale, major 
features of the dynamics can be recovered by assuming 
the pulses are ^-functions in time. With this approxima- 
tion, we are able to reduce the electron and absorption 
spectra to very simple forms in terms of the pulse areas of 
the EUV and the laser, the time delay, and the q param- 
eter, where the detailed derivations are shown in the ap- 
pendix. The modification brought by the laser will have 
maximum effects when the laser follows the EUV immedi- 
ately and its pulse area is (4iV+2)7r for integer N . In this 
case the photoelectron distribution flips about the reso- 
nance energy. For the light, the absorption line changes 
to the emission line. These dramatic changes gradually 
die down when the time delay increases in terms of the 
decay lifetime of the observed resonance. This <5-pulse 
analysis will serve to qualitatively explain the result of 
realistic numerical calculation later in this article. 



C. Propagation of light 

We have shown how to model the resonance shape of an 
AIS changed by the strong coupling laser in an atomic 



system. In principle, both the electron and absorption 
spectra carry the same essential information, depending 
on the resonance structure of the states and the laser 
parameters. However, in electron measurements, high 
precision spectroscopy is usually harder to achieve since 
the lineshape of most Fano resonances requires meV en- 
ergy resolution. In contrast, such resolution is easier with 
the transient absorption spectroscopy, as exemplified in 
a recent report (22|. Thus, in the present work we will 
put emphasis on absorption spectroscopy. 

In a linear medium, the absorption of light is commonly 
described by Beer's law as 



T(uj) = T (w)cxp [-pa(uj)L] , 



(18) 



where Tq(lu) and T{uj) are light intensities at the en- 
trance and at the exit of the medium, respectively, p is 
the number density of the particles, cr(uj) is the absorp- 
tion cross section, and L is the travel distance of light in 
the medium. This form requires that the cross section 
has no spatial and temporal dependence and is only a 
function of energy, and the transmission of light simply 
decays exponentially with its traveling distance and the 
gas density. However, in the presence of intense ultra- 
short IR pulses, the response of the medium can be more 
complicated, which should be carried out by propagating 
the exact electric field in the medium. 

In the present work, by assuming that the laser is 
loosely focused, the electric field is a function of time 
t and the spatial coordinate in the propagation direction 
z only, where the dependence on the transverse direction 
is neglected. We assume that the propagation is only 
forward in z and approximately at the speed of light in 
vacuum, c. By expressing the time in the moving frame 
t 1 = t — z/c, the Maxwell equation is reduced to 



dE(z,t') 
dz 



p dfx(z,t') 



ce 



Of 



(19) 



where p is the gas density. Because the EUV and laser 
frequencies are widely separate, they can be propagated 
separately. For the EUV, the carrier oscillation terms in 
field and in dipole are factored out, and the propagation 
for the field envelope is 



dF x {z,f) 
dz 



_P_ 

ceo 



dux{z, t') 
dt 1 



iuj x u x {z,t') 



(20) 



For the laser, the propagation simply follows Eq. (TT9")) 
where p(z,t') is replaced by ul(^') + ex.. The dipole 
oscillations corresponding to the EUV and laser fields 
are determined by Eqs. (|16j) and (|17|) . respectively, and 
these dipoles determine how the pulses evolve along z 
by Eqs. ([20]) and ([19]). The TDSE for the single-atom 
wave function and the Maxwell equation for light propa- 
gation are executed sequentially until the pulses exit the 
medium. 
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III. RESULTS AND ANALYSIS 

In the following, energy and time are expressed in elec- 
tron volts (eV) and femtoseconds (fs), respectively, and 
the signal intensities in spectra are expressed in atomic 
units (a.u.), unless otherwise specified. 

A. Results for a single atom 

In this study, the characteristics of two strongly cou- 
pled AISs are demonstrated for neutral helium, the sim- 
plest atomic system possessing electron correlation ef- 
fects. The lowest few AISs in helium are separated from 
one another on the scale of eVs and are easily coupled 
by lasers in the near IR to the visible-light energy range. 
Here we consider a laser pulse [wavelength Al — 540 nm, 
full width at half maximum (FWHM) duration tl = 9 fs, 
peak intensity II = 2 Jo where Iq = 1 TW/cm 2 hereafter 
for convenience] to resonantly couple the 2s2p( 1 P) and 
2s 2 ( 1 5) resonances, and a weak EUV SAP (central pho- 
ton energy uj x = 60 eV, FWHM duration t x = 200 as, 
peak intensity 10 10 W/cm 2 ) to excite 2s2p from the 
ground state by one-photon transition. Both pulse en- 
velopes arc in the sine-squared form. The A coupling 
scheme was chosen so that the binding energies of the 
two AISs are moderately high to prevent ionization by 
the laser. The time delay to is adjustable, which is de- 
fined by the time between the two pulse peaks, and it is 
positive when the EUV comes first. The broadband EUV 
covers roughly from 50 to 70 eV and can in principle ex- 
cite many resonances at once, but due to the choice of 
the laser wavelength and its relatively narrow bandwidth 
of about 200 meV, the coupling between the two AISs 
specified above remains our focus, where it will only be 
slightly disturbed by the presence of other states if they 
are not totally negligible. 

The time-delayed photoclcctron and photoabsorption 
spectra near the 2s2p resonance energy, calculated by 
the single-atom model, are shown in Figs.^a) andG2a), 
respectively, for the time delay ranging from —10 fs to 
40 fs. Most informative spectral features are included in 
this range. The resonant coupling laser with intensity 
2 Jo has a pulse area of 1.3 x 27T, which means that the 
laser pulse, in its whole, is able to incur Rabi oscillation 
for 1.3 cycles. 

The spectra in Figs.[2ja) and[3fa), calculated with re- 
alistic pulse shapes, retain the main features shown in 
Figs. [7] and [51 The latter two figures are derived analyt- 
ically by assuming that the two pulses are (5-functions in 
time. The most obvious differences appear when the two 
pulses overlap, where only the part of the laser after the 
EUV pulse is responsible for the Rabi oscillation. Unlike 
the (5-function short pulses, the realistic pulses bring their 
own timescales into account. Thus, we define an "effec- 
tive pulse area" by counting only the laser duration after 
the EUV peak. For to = 3.5 fs, the effective pulse area is 
2tt, i.e., the laser at that time delay supports perfectly a 
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FIG. 2: (Color online) Time-delayed photoelectron spectra 
near the 2s2p resonance excited by a 200 as EUV with a 9 
fs laser coupling, calculated by the single-atom model, (a) 
The spectra between to — —10 and 40 fs. The spectra for 
positive delays are to be compared with those in the 5-pulse 
analysis in Fig. [jj (b) The spectrum for to = 3.5 fs, which 
corresponds to the green line in panel (a). At this time delay, 
the resonance profile flips horizontally from the original Fano 
lineshape, which is shown by the gray curve. 



full cycle of Rabi oscillation. In Figs. [21(b) and[3][b), the 
spectra for to = 3.5 fs are plotted. The flipped Fano line- 
shapes discussed in Sec. Ill B 21 and in the appendix are 
seen. In particular, the electron spectrum flips horizon- 
tally where the sign of q is changed, and the absorption 
spectrum has a upside-down image, where the absorp- 
tion peak in the laser-free spectrum points downward, 
and the signals below the zero line are for the light emis- 
sion. These "flipped Fano lineshape" features have been 
shown in theoretical studies [l6l - fl8l I20I l2lj and observed 
in a recent experiment j22[. While some comprehensive 
simulations were done previously, the analytical deriva- 
tion in this work would provide a simple yet effective 
explanation that could readily be applied. 

Note that the AISs decay continuously once they are 
populated, no matter whether the laser is present or not. 
For to = 3.5 fs, at the time of the laser peak, the 2s2p 
state has already partially autoionized, and the Rabi os- 
cillation is able to affect only electrons that are not au- 
toionized yet. Furthermore, the decay of 2s 2 also "leaks" 
some electrons resulting from Rabi oscillation, by emit- 
ting photoelectrons with its own decay lifetime of 5.3 fs. 
These factors reduce the distortions that the laser would 
have imposed on the original Fano lineshape, in both the 
electron and absorption spectra. As the time delay in- 
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FIG. 3: (Color online) Same as Fig. [2] but for the photoab- 
sorption spectra. The spectra for positive delays are to be 
compared with the spectra in the 5-pulse analysis in Fig. [8] 
The resonance profile for to = 3.5 fs flips vertically from the 
original Fano lineshape shown by the gray curve. Note that 
the negative signals at to = 3.5 fs around the resonance energy 
represent the emission of light. In contrast, the photoelectron 
spectra in Fig. [2] do not have negative signals. 



FIG. 4: (Color online) Time-delayed photoabsorption spectra 
of a 200 as EUV with a 50 fs dressing laser in the single-atom 
model, (a) The spectra between to — —50 fs and 80 fs. The 
EIT condition is controlled by the time delay between the 
two pulses, (b) The spectrum for to = 10 fs shows the largest 
separation in the Autler-Townes doublet, which matches the 
Rabi frequency at the laser peak. 



creases, more electrons are autoionized before the laser, 
and less are coupled by the laser to participate in the 
Rabi oscillation, and the resonance profile gradually ap- 
proaches the original Fano lineshape. 

While similar three-level coupling schemes have been 
used for the electromagnctically induced transparency 
(EIT) effect [H,^, the spectral features in a typical EIT 
setup, such as the transparency line or Autler-Townes 
doublet [13], are not recovered by our photoabsorption 
calculation in Fig. [3] Evidently, the approximation with 
(5-function pulses breaks down in the long dressing field. 
To differentiate our mechanism from EIT, we make an 
additional calculation by keeping all the parameters, but 
the laser duration is extended from 9 fs to 50 fs. The 
resultant spectra are plotted in Fig. [4] It shows that the 
long laser pulse splits the resonance peak into two peaks, 
with the maximum separation 0.45 eV at to = 10 fs. The 
separation is exactly the Rabi frequency at the laser peak, 
where the Autler-Townes doublet is reproduced. As the 
time delay changes from this optimal value, the separa- 
tion decreases, and finally drops to for to < —40 fs 
or to > 70 fs. It shows that contrary to short dress- 
ing pulses, long dressing pulses with many Rabi cycles 
creates the condition satisfying the "dressed state" pic- 
ture behind the EIT phenomenon, and the split states 
appear as expected. By controlling the overlap between 



the pulses, the doublet can be "turned on" or "turned 
off." This is the theme of some recent studies applying 
dressing lasers of tens to hundreds of femtoseconds in 
an EIT scheme while measuring the time-delayed EUV 
transmission (l3l. [TBI. l28[j . 

B. Results for a gaseous medium 

In Sec. 1111 Al the photoabsorption spectrum reveals 
the appearance of emission line at the resonance for 
to = 3.5 fs. However, it is physically impossible for the 
light at certain frequency to gain intensity indefinitely 
when passing through a medium. A realistic pulse has 
no singularity in the frequency domain. Thus, by intu- 
ition, the single-atom result does not represent what will 
actually be measured in a dense gaseous medium. With 
the same EUV and laser pulses used in Sec. lIII Al we con- 
sider a gaseous medium made of non-interacting helium 
atoms with number density p = 8 x 10 17 cm~ 3 (equiv- 
alent to pressure of 25 Torr at room temperature) and 
thickness L = 2 mm, and calculate the transmitted light 
spectra at the exit of the medium. 

The transmitted EUV is plotted in Fig.[5]for three laser 
intensities, II = I.IIq, 4.5/oi an d 10/o- where II is the 
peak intensity. They safely fall into the intensity range 
considered by our model. The EUV spectra with overlap- 
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ping pulses fixed at to = are plotted in Fig. [5ja) . The 
effective pulse areas for the three intensities are n, 2n, 
and 3-7T, which are responsible for the half, one, and one- 
and-half cycles of Rabi oscillation. For l.lirj) the spec- 
trum is mostly flat, indicating that while the laser moves 
the electron population from 2s2p to 2s 2 , the 2s2p AIS 
"disappears," leaving only the directly ionized photoelec- 
trons in the vicinity without the hint of autoionization. 
For 4.5/o, the laser drives the electrons back to 2s2p with 
a phase shift of tt. The emission line that has been dis- 
cussed in the single-atom picture is observed here, where 
the spectral peak appears to be more than 50% of the 
signal intensity of the incident pulse, which is effectively 
a partial enhancement of the SAP passing through the 
gas. When the laser intensifies further to 10 Io, the Rabi 
oscillation drives the electrons away from 2s2p again; the 
enhancement is gone, and the spectrum looks flat other 
than a slight bump. These three laser intensities demon- 
strate how the number of cycles in the Rabi oscillation 
controls the phase of the 2s2p AIS, which then determines 
the resonance lincshape at the end of the autoionization. 

The presentation in Fig. [5ja) is for the overlapping 
pulses, i.e., the Rabi oscillation happens within the laser 
duration at the beginning of the 17 fs decay of 2s2p. 
It is a simplified case where autoionization occurs ap- 
proximately after the population transfer by the laser is 
done. However, for larger time delays, the laser comes 
after 2s2p decays for some time and only affects the later 
part of the decay. The total wave packet is the coherent 
sum of autoionization before the laser, and the quantum 
path going through the Rabi oscillation. In Fig. [5jb), 
the total signal yields in the 50 meV range around the 
2s2p resonance for different intensities and time delays 
are plotted. This plot is equivalent to a measurement 
of the EUV signal at the resonance energy with a spec- 
trometer of 50 meV resolution in order to investigate the 
enhancement. For each intensity, the yield oscillates with 
the time delay. As seen in the figure, 4.5io is n0 longer 
the definite optimal intensity to enhance the transmitted 
SAP resonance. For example, when the 10/o pulse is at 
to = —2 fs, the enhancement is higher than what can be 
achieved by the 4.57o pulse. This is because even when 
both intensities have the same effective pulse area, the 
lO/o pulse accomplishes the Rabi oscillation in a slightly 
earlier stage in the decay of 2s2p, where less electrons 
escape its influence by autoionization. It is very remark- 
able that this partial enhancement of a SAP is controlled 
within an energy scale of about 100 meV in the form of 
resonance lineshape, and within a timescale of about 1 fs 
in the form of time delay. 

It has been demonstrated in Fig. Ufa) that while the 
resonance lineshape of the EUV pulse is controlled by 
the laser coupling, the spectral profile away from the res- 
onance energy always drops in light transmission spec- 
troscopy. It is imaginable that the overall light attenuates 
when propagating in the gas. It is thus very intriguing to 
know how the exceptional enhanced resonance peak holds 
its shape and strength in the propagation. In Fig. [BJ we 
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FIG. 5: (Color online) Transmission of EUV through a helium 
gaseous medium dressed by a 9 fs laser of peak intensities (7l) 
of (laser free), l.Uo, 4.57~o, and 10/o- The incident light is 
plotted as a reference, (a) Transmitted EUV profiles for over- 
lapping pulses (to = 0), where the effective pulse areas of the 
laser are 0, 7T, 2tt, and 37r. The flipping of the resonance pro- 
file is predicted as the main feature by the single-atom case, 
although more structures are seen here. The signal at the 
resonance for II = 4.57o is enhanced significantly, with more 
than 50% increase from the incident light. This measurable 
quantity can only be obtained by applying field propagation, 
(b) Total photon signals (yield) gathered in the 50-meV en- 
ergy range at the resonance [indicated in panel (a) by the 
yellow vertical stripe] as functions of the time delay. This en- 
ergy range complies with the spectrometer resolution under 
current technology. The yield for each laser intensity oscillates 
with the time delay, which shows a high degree of control in 
the time domain. 



show the transmitted EUV spectra at different distances 
along the beam path with the 9-fs, 4.5/o dressing pulse at 
time delay. The transmitted EUV intensities and the 
corresponding effective absorption cross sections, defined 
by [sec Eq. (Tg])] 



(T off (w) = -In — — - 

pL T {u) 



(21) 



are shown in Figs.UJa) and (b), respectively. With such 
laser parameters, the "wings" on the two sides of the 
resonance profile descend, but the enhancement peak is 
up to about 50% of the incident intensity, and roughly 
maintains its height from 1 to 3 mm. The persistence of 
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the enhancement peak, apart from the exponential decay 
of the rest of the signals along the propagation, displays 
an impressive nonlinear response of the gas to the laser 
pulse. This violation to Beer's law is in a very selective 
energy range and is sensitive to the parameters of the 
dressing laser. To elucidate the nonlinearity of the dress- 
ing laser, <j c r{ui) without the dressing field is plotted in 
Fig. [SJc). In this case, the propagation of the EUV ac- 
tually follows Beer's law stated in Eq. (fT8|) . where the 
cross section determined by the original Fano lineshape 
remains the same over the distance. The pulse over the 
whole energy range is absorbed indifferently 



IV. CONCLUSION 

A three-level autoionizing system coupled by a SAP 
and a time-delayed intense femtosecond laser have been 
modeled. The photoclcctron and photoabsorption spec- 
tra are responsive most obviously to the laser intensity 
and its time delay. We analyze the spectra by assuming 
both pulses are infinitely short, where the responses can 
be explained by the precise control of the Rabi oscillation 
between the two AISs. The analysis also suggests that 
the phase change in the wave packet during the Rabi 
oscillation is mapped out in the flipping of the Fano line- 
shape. The model calculation with realistic pulse pa- 
rameters also produces these significant features. In con- 
trast, with long dressing fields, these spectral features 
are replaced by the Autler-Townes doublet, which is well 
studied in systems coupled by stable lasers. This contrast 
demonstrates that the ultrashort coupling demolishes the 
dressed-state picture, and the dynamics of the system is 
better viewed directly in the time domain. 

Focusing on the response of the SAP light in this 
coupled system, we further incorporate the single-atom 
response with the pulse propagation in the gaseous 
medium. The main controllable features of the resonance 
profile shown by a single atom are well preserved in the 
gaseous medium. By applying a laser pulse with a 2tt ef- 
fective pulse area to couple the 2s2p and 2s 2 resonances 
in the dense helium gas, an enhanced peak at the 2s2p 
resonance shows up in the transmitted SAP spectrum, 
where the signal intensity retains more than 50% gain 
from the entrance to the exit of the medium. At the 
same time, the signals away from the resonance energy 
drops along the propagation as described by Beer's law. 
It can be viewed as the pulse shaping of a SAP that is 
measurable in the energy domain. This result illustrates 
a strong nonlinear manipulation of an SAP by an ultra- 
fast dressing field, which would open new possibilities of 
control of attosecond dynamics. 
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FIG. 6: (Color online) (a) EUV transmission profiles for prop- 
agation lengths L = 1, 2, and 3 mm in the gaseous medium 
with the overlapping 9 fs, 4.57o dressing laser pulse. When 
the dressing laser is turned on, the enhanced peak persists 
while the surrounding background attenuates. Panels (b) 
and (c) show effective absorption cross sections a e ff(cj) for 
the same propagation lengths as in panel (a) with and with- 
out the dressing laser, respectively. With the dressing laser, 
<j e ff(ui) changes with L at the resonance but stays the same 
in the outside region. Without the dressing field, a e g(uj) over 
the whole energy range does not change with L, where Beer's 
law applies. The contrast between the cases with and without 
the dressing laser demonstrates a strong nonlinear response of 
the gas to the laser field which is reflected in the transmission 
of the attosecond pulse. 
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In the short-pulse limit, we assume that the field en- 
velopes are rectangular with given carrier frequencies. As 
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the durations approach zero, the envelopes are 
Fx{t) = ^D gbl S(t) 
F L (t) = ^D blb2 S(t-t ), 



(Al) 



where Ax and Al arc equivalent to the pulse areas for 
the EUV and laser, respectively, and to is the time delay 
between the pulse peaks, which is positive when the EUV 
comes first. Since the intensity of the laser is not enough 
to excite the atom from the ground state, we only con- 
sider the time delay to > 0. The infinitely short pulses 
cut the time into regions where the TDSE in each region 
can be solved analytically. The EUV is weak so that 
perturbation is applied on its transition. In order to ob- 
serve the universal features of the electron and absorption 
spectra regardless of the atomic parameters, the energy 
is scaled by e = 2(E — E bl ) /Ti, where Ti is the resonance 
width (see Sec. [TTJ) . Note that now E bl represents either 
the electron energy or the photon energy depending on 
whether we look at the electron or absorption spectra. 
This scaling is the same as Fano's characterization of 



resonances 



24j . The time delay is scaled by r = Tito/2. 



With some algebra, the electron spectrum is found to 
be given by 



P(e) 



1-4 



A' 



where 



and 



[( q + e) 2 -ce- T f P (e)], (A2) 



1 - cosOU/2) (A3) 



f P (e) = 2(q + e) [q cos(er) + sin(er)] +ce' T (q 2 + 1). (A4) 
Similarly, the absorption spectrum is given by 



S(e) = ^[(q + er-ce- T f s (e)} 



(A5) 



where 



fs(e) = (q 2 -l + 2qe) cos(er) + [2q + (1 - q 2 )e] sin(er). 

(A6) 

Both spectra are composed of two terms. The first term 
is the original Fano lineshape defined by the \b\)-\Ei) res- 
onance alone, without any influence from the laser. The 
second term decays exponentially with the time delay t, 
and is proportional to the coefficient c determined by the 
laser pulse area Al as defined by Eq. (|A3[) . The minimum 
value of c is 0, which corresponds to Al = 0,4ir, 8ir, 
where the Rabi oscillation has even number of cycles. In 
this condition, the electrons coming back to carry no 
additional phase, and the laser effectively does nothing. 
On the contrary, the maximum value of c, 2, corresponds 
to Al = 2-7T, 6-7T, 107r, ... and odd number of cycles. The 
returning electrons have additional phase of ir. 

In order to demonstrate the influence of laser strength 
more clearly, the electron and absorption time-delayed 



spectra for Al = 2tt are plotted in Figs. [7] and [H re- 
spectively As seen in both figures, when r is large, very 
little change is brought by the laser onto the original Fano 
lineshape. However, in the r — > limit, the electron spec- 
trum flips in the energy dimension around e = 0, and the 
absorption spectrum displays a strong emission line near 
e = 0. 



CD 
T3 



FIG. 7: Photoelectron spectra with infinitely short pulses, as 
described by Eq. (| A2 j> ■ The energy is scaled by Fi/2, and 
the time delay is scaled by 2/Ti. The pulse area of laser is 
Al = 2tt, which indicates a full cycle of Rabi oscillation by 
the laser coupling. When the time delay is large, the spectrum 
approaches the original Fano lineshape. When the time delay 
is 0, the spectrum is close to a horizontally flipped image of 
the original Fano lineshape. 



The significant change of spectral features at r = 
can be expressed analytically by reducing Eqs. (|A2[) and 
[!) to 



P(e) = \A x [ 
S(e) = \Ax\- 



(-g + e) 2 +4 
1 + e 2 ' 



1 



(A7) 
(A8) 



In the electron spectrum, the mirror image of Fano res- 
onance, indicated by the minus sign of the q parame- 
ter, is added onto a Lorentzian shape. In the absorption 
spectrum, beside an additional constant background, the 
whole resonance function flips upside-down, which means 
that the absorption line becomes the emission line. Note 
that Eqs. (fA7|) and (|A"8) arc the results for a full Rabi 
cycle and represent the extreme changes that the strong 
coupling can bring. For non-integer number of cycles, 
among the two resonances, some electrons will stay at 
\b 2 ) at the end of the laser, and will not be combined 
with the photoelectrons in the neighborhood of \bi). 
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FIG. 8: Same as Fig. [7] but for the photoabsorption, as de- 
scribed by Eq. (|A5I) . When the time delay is 0, the spectrum 
is close to a vertically flipped image of the original Fano line- 
shape. 
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